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ABSTRACT 

Recent observations have revealed that the time evolution of the SO number frac- 
tion at intermediate and high redshifts (0.2 < z < 0.8) is more dramatic in groups 
of galaxies than in clusters. In order to understand the origin of SOs in groups, we 
investigate numerically the morphological transformation of spirals into SOs through 
group-related physical processes. Our chemodynamical simulations show that spirals 
in group environments can be strongly influenced by repetitive slow encounters with 
group member galaxies so that those with thin disks and prominent spiral arm struc- 
tures can be transformed into SOs with thick disks and without prominent spiral 
arm structure. Such tidal interactions can also trigger repetitive starbursts within the 
bulges of spirals and consequently increase significantly the masses of their bulges. 
Owing to rapid consumption of gas initially in spirals during the bulge growth, the 
SOs can become gas-poor. The SOs transformed from spirals in this way have young 
and metal-rich stellar populations in the inner regions of their bulges. The simulated 
SOs have lower maximum rotational velocities and flatter radial line-of-sight velocity 
dispersion profiles in comparison to their progenitor spirals. The formation processes 
of SOs due to tidal interactions depend not only on the masses and orbits of the pro- 
genitor spirals, but also on group mass. A significant fraction (10 — 30%) of stars and 
gas can be stripped during this spiral to SO morphological transformation so that in- 
tragroup stars and gas can be formed. Based on these results, we discuss structures, 
kinematics, chemical properties, and the Tully-Fisher relation of SOs in groups. 

Key words: galaxies: evolution- galaxies: elliptical and lenticular, cD- galaxies: 
high-redshift- galaxies: starburst- galaxies: bulges 



1 INTRODUCTION 

Many authors have long discussed when and how SOs were 
formed in different environments based on their observed 
fundamental properties such as apparent shapes (e.g., van 
den Bergh 1976; Holden et al. 2009), gas content (e.g., 
Bothun 1982; Welch & Sage 2003), stellar structure and 
kinematics (e.g., Kormendy & lUingworth 1982; Burstein 
1979; Fisher 1997; Emsellem et al. 2007; Laurikainen et 
al. 2010), counter-rotating components (e.g., Bertola et al. 
1992), stellar population ages (e.g., Kuntschner & Davies 
1998), color-magnitude relation (e.g., Ellis et al. 1997; 
Stanford et al. 1998; van Dokkum et al. 1998), and the 
morphology-density relation (e.g.. Dressier 1980). Most of 
these studies considered that SOs originate from spirals and 
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thereby discussed what physical roles galaxy environments 
(e.g., groups and clusters) play in transforming spirals into 
SOs. Recent observations have suggested that not only en- 
vironment but also galaxy mass might be important in SO 
formation (e.g., Tasca et al. 2009; Vulcani et al. 2010). The 
origin of SOs has also been linked observationally to the evo- 
lution of other galaxy populations such as dusty starburst 
galaxies (e.g., Geach et al. 2009) and poststarburst ones with 
"E+A" spectra (e.g., Pracy et al. 2009). 

A growing number of observational studies have re- 
vealed a smaller fraction of SO galaxies in distant clusters of 
galaxies (relative to that observed in present day clusters), 
suggesting dramatic morphological evolution of spirals into 
SOs within these environments (Dressier et al. 1997; Couch et 
al. 1998; van Dokkum et al. 1998; Fasano et al. 2000; Smith 
et al. 2005; Postman et al. 2005; Desai et al. 2007; Poggianti 
et al. 2008; Vulcani et al. 2010), though such SO evolution 
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may not be real but due to morphological classification er- 
rors (e.g., Andreon 1998). Recent photometric studies of the 
structure of SOs have revealed intriguing scaling-relations be- 
tween their bulges and disks and differences in bar strength 
between spirals and SOs (e.g., Buta ct al. 2010; Laurikainen 
et al. 2010). These observations have provided new clues £is 
to whether and how SOs are formed from spirals in different 
environments. 

A number of theoretical models for the transforma- 
tion of spirals into SOs have been proposed, which includes 
ram pressure stripping (e.g., Gunn & Gott 1972; Farouki & 
Shapiro 1980; Abadi et al. 1999), tidal encounters (e.g., Icke 
1985), tidal compression by the gravitational field of clusters 
(e.g., Byrd & Valtonen 1990), truncation of gas replenish- 
ment (e.g., Laxson et al. 1980; Bekki et al. 2002), and minor 
and unequal-mass merging (e.g., Bekki 1998). Given that 
SOs exist in different environments such as the field, groups, 
and clusters, with the SO fraction depending on local galaxy 
density (e.g.. Dressier 1980), all these physical mechanisms 
could be important and their relative importance of would 
depend on environmental parameters such as the total group 
and cluster mass. Thus it remains unclear which theoretical 
model is the most important in SO formation for a given 
environment. 

Recent observational studies on the evolution of the SO 
fraction in groups and clusters have found that SO evolu- 
tion is significantly more dramatic in groups than in clusters 
(e.g., Wilman ct al. 2009; Just et al. 2010). These observa- 
tions suggest that cluster-related physical processes such as 
ram pressure stripping (e.g., Gunn & Gott 1972) are not 
responsible for the formation of the majority of SOs. Also, 
it is well known that galaxy bulges in SOs arc systemat- 
ically more luminous than those in spirals (e.g., Dressier 
1980; Simien & de Vaucouleurs 1986; Christlein & Zablud- 
off 2004), which implies that bulges in spirals need to grow 
significantly in order to bo transformed into SOs: disk-fading 
alone, due to the truncation of star formation, cannot be the 
main mechanism of SO formation. Thus galaxy interactions 
and merging, which axe highly likely to occur in groups, are 
a promising mechanism for SO formation and thus need to 
be investigated theoretically using numerical simulations. 

The purpose of this paper is thus to investigate how SOs 
are formed in the group environment using chemodynami- 
cal numerical simulations that enable us to investigate both 
dynamical properties (e.g., stellar kinematics) and star for- 
mation and chemical evolution histories for spirals within 
groups in a fully self-consistent manner. We consider that 
both (i) slow tidal encounters of galaxies with relative veloc- 
ities (K) less than 400 — 500 km s~^, and (ii) group tides arc 
responsible for the transformation of spirals into SOs in the 
group environment. We particularly investigate how spirals 
with smaller bulges are transformed into SOs with bigger 
ones as a result of the tidal fields within groups and group 
member galaxy interactions. We discuss how this transfor- 
mation process depends on the masses and orbits of galaxies, 
the interaction histories of group members, and the physical 
properties of groups such as their masses and sizes, based 
on a large and systematic parameter study of SO formation. 

The plan of the paper is as follows: In the next section 
we describe our numerical model for SO formation in groups. 
In §3, we present the numerical results mainly on the phys- 
ical properties of the simulated SOs. In §4, we discuss these 



results in the context of key observational results on the 
physical properties of SOs. We summarize our conclusions in 
§5. 



2 THE MODEL 

2.1 Group of galaxies 

In order to simulate the time evolution of spirals in groups, 
we use the latest version of GRAPE (GRavity PipE, 
GRAPE- 7), which is the special-purpose computer for grav- 
itational dynamics (Sugimoto et al. 1990). We use our 
original GRAPE-SPH code (Bekki 2009) which combines 
the method of smoothed particle hydrodynamics (SPH) 
with GRAPE for calculations of three-dimensional self- 
gravitating fluids in astrophysics. In the present models, 
both a spiral galaxy and the dark matter halo of its host 
group of galaxies are represented by N-body particles so 
that not only dynamical friction of the spiral against the 
dark matter halo of the group but also dynamical influences 
of the halo on the spiral can be self-consistently investigated. 
Group member galaxies other than the spiral arc modeled 
as point-mass particles, because we consider that such mod- 
eling is enough to grasp the essential ingredients of tidal 
influences of group member galaxies on the evolution of the 
spiral in the present study. 

The structure of a group is modeled using an "NFW" 
profile predicted by the cold dark matter cosmology 
(Navarro et al. 1996) as follows: 



p(r) : 



(r/r,)(l + r/r.)2' 



(1) 



where r, po, and rs are the spherical radius, the characteristic 
density of a dark halo, and the scale length of the halo, 
respectively. The mass and size of a group are represented 
by Mgr and Rgr, respectively. The c parameter (= Vs/rviv, 
where rvir is the virial radius of the NFW profile) for a 
group with Mgr is chosen according to the predicted c-Mgr 
relation in the ACDM simulations (e.g., Neto et al. 2007) 
:here the dark matter mass of the group corresponds to Mgr. 
A reasonable value of c is thus 5.6 for Mgr = 2 x lO^^M© and 
4.7 for Mgr = 10^*Mq. We mainly show the results of the 
models with Mgr = 2 x lO^^M© and i?gr = 535 kpc, though 
we investigate different models with Mgr ranging from 5 x 
IO^^Mq to 10"Mq. 

Galaxies in a group are represented by collisionless par- 
ticles and their spatial distribution follows the NFW profile 
with c = 3, which is consistent with recent i^-band observa- 
tional studies on galaxies distributions in groups and clus- 
ters (e.g., Lin et al. 2004). The canonical Schechter function 
(with slope of —1.07) is adopted for generating a galaxy lu- 
minosity/mass function for luminosities ranging from O.OIL* 
to 2.5//* in a group. We assume that the mass-to-light-ratio 
(here mass includes dark matter) for each individual galaxy 
is 20 and thereby calculate the total mass according to the 
allocated luminosity. The total mass (thus rmmber) of galax- 
ies in a group with Mgr is determined by the mass-to-light- 
ratio, that itself is dependent on Mgr (Maxinoni & Hudson 
2002); 



Mgr/Lgr = 350( 



AL 



5 X 1O"M0 



(2) 
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Table 1. Description of the fiducial model and the range of model parameters for other representative models. 

Model Afd^ fidl /b^ fg£ ^gr^ Orbital types /_ i?ini£, /vj^, by j_ 

Fiducial 6 X 17." O.lT 0." 2 X l"^ Normal ~ 4rs~ 0.5~ -~ 

Representative (1.2-60) xlO^ 6.6-17.5 0-0.7 0-0.5 (0.05-1.0) X lO" Normal/Infall (1-7)xt-s 0.2-1 (0.5-2) xr^ 



The initial disk mass of a spiral in units of Mq. 
^ The initial size of a stellar disk for a spiral in units of kpc. 
^ The initial bulge mass fraction (M^/Md) of a spiral. 

The initial gas mass fraction (Mg/M^) of a spiral. 
^ The initial total mass of a group of galaxies in units of Mq . 

f "Normal model" means that a spiral can orbit the center of its host group whereas "Infall model" means that a spiral can infall onto 
a group from the outside of the group. 

9 The initial distance of a spiral from the center of its host group. This iJinj is the x-position of the spiral (and alwrays positive) and 

represents the scale-length of the adopted NFW dark matter halo of the group. 
^ For a normal model, the initial i;-component velocity of a spiral is given as fvVc, where Vc is the circular velocity at the position of 
the spiral in a group. For an infall model, the x-componcnt velocity of a spiral is given as —fvVc (because Ri^i is always positive in the 
present study). 

* The impact parameter of the orbit of a spiral in a infall model. This by corresponds to the initial i;-component of the location of the 
spiral. 



where Lgr is the total luminosity of galaxies in the group. 
Therefore, firstly Lgr for a group with mass Mgr is derived 
from the above equation (2) and then the total number of 
galaxies for the group is determined from the derived Lgr 
and the adopted luminosity function of galaxies. As an ex- 
ample, a group with Mgr — 2 x lO^^M© has 87 galaxies. 
Group member galaxies in a group have an isotropic veloc- 
ity dispersion just as the dark matter of the cluster does. 



2.2 Spiral galaxy 

The total mass and the size of the disk of a spiral galaxy 
with total galaxy mass Mgai (inclusive of dark matter) are 
Md and Rd, respectively. Henceforth, all masses and lengths 
are measured in units of Md and Rd, respectively, unless 
specified. Velocity and time are measured in units of n = 
{GMd/Rdy^'^ and tdyn = (Rl/GMdY^'^, respectively, where 
G is the gravitational constant and assumed to be 1.0 in 
the present study. If we adopt Aid ~ 6.0 x 10^° M© and 
Rd = 17.5 kpc as a fiducial value, then v — 1.21 x 10*^ 
km/s and tdyn = 1-41 x 10* yr, respectively. The disk is 
composed of a dark matter halo, a stellar disk, a stellar 
bulge, and a gaseous disk. Gaseous halos that were included 
in our previous works (e.g., Bekki 2009) are not included in 
the present paper, because we do not discuss ram pressure 
stripping of halo gas within groups. 

The mass ratio of the dark matter halo to the stel- 
lar disk in a spiral model is fixed at 9 for most mod- 
els and the density distribution of the halo is represented 
by the NFW profile. We consider that a reasonable value 
of the c-parameter is 7.8 (consistent with the results by 
Neto et al. 2007) for Milky Way-type disk galaxies with 
Md ~ 6 X 1O^°M0. The value of rg is chosen such that the 
rotation curve of a disk is reasonably consistent with obser- 
vations: the majcimum circular velocity of a Milky Way-like 
disk galaxy in the fiducial model is 245 km s~^. 

The bulge of a spiral has a mass of Afb, a size of 
Rh and a scale-length of i?b,o and is represented by the 
Hernquist density profile. The bulge is assumed to have 
isotropic velocity dispersion and the radial velocity disper- 
sion is given according to the Jeans equation for a spherical 



system. The bulge-mass fraction (/b = Mb/Md) is a free 
parameter. We mainly investigate "Milky Way" models (re- 
ferred to as "MW" from now on) in which /b — 0.167 and 
Rh = 0.2i?d (i-e., Rh.o = 0.047?d)- In order to determine 
Rh for a given Mb and fh, we use the Faber- Jackson rela- 
tion (L oc (J*, where L is the total luminosity of a galaxy 
and a is the central velocity dispersion; Faber & Jackson 
1976) and the virial theorem. The "B/T ratio" is defined 
as B/T — Mh/{Md + Mb), which means that it is different 
from fh- The initial and final B/T are referred to as {B/T)i 
and {B/T)i, respectively. 

The radial (R) and vertical (Z) density profiles of the 
stellar disk are assumed to be proportional to exp(— i?/_Ro) 
with scale length Ro — 0.2 and to sech^ (Z/Zq) with scale 
length Zo — 0.04 in our units, respectively. In addition to 
the rotational velocity caused by the gravitational field of 
disk, bulge, and dark halo components, the initial radial 
and azimuthal velocity dispersions are assigned to the disc 
component according to the epicyclic theory with Toomre's 
parameter Q = 1.5. The vertical velocity dispersion at a 
given radius is set to be 0.5 times as large as the radial 
velocity dispersion at that point, as is consistent with the 
observed trend of the Milky Way (e.g., Wielen 1977). 

We investigate models with different Mgai and adopt 
Freeman's law (Freeman 1970) to determine Ro of a disk 
galaxy according to its disk mass: 



Ro = Cd( 



Md 



6 X IO^oMq 



0.5 

) kpc. 



(3) 



where Cd is a normalization constant. We consider that 
Cd ~ 3.5 kpc is a reasonable value, because it is consis- 
tent with the disk structure of the Galaxy. We thus adopt 
Cd = 3.5 kpc as a standard value for luminous disk galax- 
ies. However, as shown in Kauffmann et al. (2003), low- 
luminosity disk galaxies have low surface stellar densities 
and thus Rq determined by the above equation would not be 
so appropriate (i.e., significantly smaller than the observed). 
We thus adopt Cd ~ 6.2 kpc and 9.3 kpc to model low- 
luminosity and low-surface-brightness galaxies. The galaxies 
with Cd = 3.5 kpc are referred to as "high-surface-brightness 
galaxies (HSBs)" and those with Cd = 9.3 kpc are as "low- 
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surface-brightness galaxies (LSBs)" just for convenience. For 
example LSBs with Md = 1.2 x lO^'M© have Ro = 1.3 kpc for 
Cd = 9.3 kpc. Structural and kinematical properties of daxk 
matter halos and stellar disks are assumed to be self-similar 
between models with different Mgai. 

2.3 Gas, star formation, and chemical evolution 

Following our previous models for galactic global star for- 
mation from gas in galaxies influenced both by galaxy in- 
teractions and tidal fields of their host larger halos (Bekki 
& Chiba 2005) and chemical evolution and supernova feed- 
back effects (Bekki & Sliioya 1999), wc model physical pro- 
cesses of star formation and chemical evolution as follows: 
The gas mass fraction (/g) is assumed to be a free parameter 
ranging from 0.1 to 0.5. An isothermal equation of state is 
used for the gas with temperatures of 10* K for models with 
Md = 6 X 10^°Mq and the initial temperature (Tiso) of disk 
gas is assumed to depend on Mgai. The gas disk of a spiral is 
represented by SPH particles and has an exponential radial 
density profile that is exactly the same as the stellar disk 
has. 

In order to construct as realistic a gas disk model as 
possible, we consider the radial dependence of the gas mass 
fraction Fg{r) in the initial disk. We expect that the inner 
gas mass fraction in the disk is smaller than the outer value 
owing to more rapid consumption of gas in the inner regions 
with higher gas density. We therefore adopt the following 
rule: 

Fg(r)cxW(r)<x5£M<xEg'^(r), (4) 
Sg(r) 

where r, tsf, Sg, Eg, and a are the distance from the center 
of the disk, the gas consumption time scale, the initial gas 
density, the gas consumption rate, and the parameter con- 
trolling the radial dependence. Since we adopt the Schmidt 
law (Schmidt 1959) with exponent of 1.5 for star formation 
(described below), a reasonable value of a is —0.5. 

According to the value of Fg{r) at each radius derived 
from the above equation (4), we determine a reasonable 
number of stellar and gaseous particles at each radius and 
thereby allocate these particles to each radial bin. By as- 
suming that the disk is composed only of gas initially, we 
determine Fg{r) through equation (4) and derive a reason- 
able radial distribution of gas and stars. If we adopt models 
with constant -fg(r), then a strong initial starburst can oc- 
cur in the central region of a gas-rich spiral. We therefore 
consider that this starburst in an isolated spiral is not rea- 
sonable and realistic. Furthermore the adopted models with 
radial dependence in Fg (with no initial starbursts) enable us 
to clearly understand how galaxy interactions trigger star- 
bursts within the bulges of spirals. 

Star formation is modeled by converting the coUisional 
gas particles into coUisionless new stellar particles accord- 
ing to the algorithm of star formation described below. We 
adopt the Schmidt law with exponent 7 = 1.5 (1.0 < 7 < 
2.0, Kennicutt 1998) as the controlling parameter of the rate 
of star formation. The amount of gas consumed by star for- 
mation for each gas particle in each time step is given as: 

Pg oc pg^, (5) 

where pg is the gas density around each gas particle. These 



stars formed from gas are called "new stars" (or "young 
stars") whereas stars initially within a disk are called "old 
stars" throughout this paper. 

Chemical enrichment through star formation and super- 
nova feedback during the evolution of the disk is assumed 
to proceed both locally and instantaneously in the present 
study. We assign the metallicity of the original gas parti- 
cle to the new stellar particle and increase the metallicity 
of each neighboring gas particle with the total number of 
neighbor gas particles equal to A^gas, according to the fol- 
lowing equation that governs chemical enrichment: 

AMz = {^>i?metms + (1.0-ftnot)(1.0-^0msymet}/iVgaB(6) 

where AAfz represents the increase in metallicity for each 
gas particle. Zi, Rmct, m^, and j/mot in the above equation 
represent the metallicity of the new stellar particle (or that 
of original gas particle), the fraction of gas returned to the 
interstellar medium, the mass of the new star, and the chem- 
ical yield, respectively. 

Since our main focus is morphological transformation 
of luminous spirals (like the Galaxy) into SOs, we mainly 
investigate chemical evolution of the MW models. For mod- 
els with Md = 6 X lO^^Mg, the values of Rmct, J/mot, and 
initial metallicity are set to be 0.3 and 0.01, and 0.01 (i.e., 
[Fe/H]=— 0.3), respectively. For comparison, we also inves- 
tigate models with the initial metallicities of 0.02 (the solar 
abundance). We adopt the same model for supernova feed- 
back effects (on ISM of galaxies) as used in our previous 
chemodynamical simulations (Bekki & Shioya 1999). In the 
present study, all of the energy from a supernova (10®^ erg 
per a supernova) is converted into kinematic energy of gas 
around the supernova. 

2.4 Orbits of spirals within groups 

We investigate the dynamical and chemical evolution and 
star formation history of a spiral galaxy in a group for 
a given model: only one galaxy is modeled fully self- 
consistently by N-body particles and other galaxies orbit- 
ing the group are represented by point-mass particles with 
the spherical Plummer softening lengths. The orbit of our 
model spiral is assumed to be influenced both by the gravi- 
tational potential resulting from the dark halo component of 
the group and by the group member galaxies. The adopted 
group potential is spherically symmetric (not triaxial), and 
thus the initial orbital plane of the spiral is set to be the x-y 
plane (=disk plane of the spiral) in all our models. Owing 
to dynamical friction of the spiral against the dark matter 
halo of the group, the spiral can sink into the central region 
of the group if the total mass of the spiral is large enough. 

The center of the group is always set to be {x,y,z) = 
(0,0,0) whereas the initial position of the spiral is set to 
be {x,y,z) = (J?ini, 0, 0). The initial velocity of the spiral 
{vx,Vy,Vz) is set to bo (0, fvVc, 0), where /v and 14 are the 
parameters controlling the orbital eccentricity (i.e, the larger 
/v is, the more circular the orbit becomes) and the circular 
velocity of the group at R = iiini, respectively. For these 
models, Vs < Rini ^ Tvir and 0.25 ^ /v ^ 1.0, and thus 
spirals can initially orbit the centers of their host groups. 
These models with J?ini ^ rvir are referred to as "normal 
models" . 

We also investigate "infall models" in which spirals are 
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initially located outside the virial radius of their host groups. 
For these models, i?ini is larger than rvir (and y is not 0) 
and Wx = fvVc (i.e., Vx is not and instead Vy = 0). In 
infall models, spirals move initially parallel to the x-axis 
toward the negative x direction and the initial value of y 
corresponds to the impact parameter (6y) of the orbit. For all 
infall models, i?ini is set to be 7rs and the impact parameter 
by and fy are the two key parameter that determine the 
orbits within groups. 

2.5 Choice of parameters 

2. 5. 1 Fiducial model 

We run models with different Md, -Rd, /b, /g, M^r, Rini, 
and /v in order to investigate the formation processes of SOs 
and their dependences on model parameters. The range of 
each model parameter is shown in the Table 1. We mainly 
show the results of the "fiducial model" in which Ma = 6.0 
X 10^" Mo i?d = 17.5 kpc (i.e., Cd = 3.5 kpc), /b = 0.17 
(i.e., MW model), /g = 0.1, Mg^ = 2 x IO^^Mq, ilini = An, 
and /v = 0.5 (i.e., "normal model"). This is mainly because 
this model shows the typical behavior of SO formation for 
spirals in groups. We run a model for ~ 6 Gyr in order to 
investigate whether the spiral can be transformed into an SO 
in the model. 

2.5.2 Multiple gravitational softening lengths 

Since we need to run ~ 100 models for allocated compu- 
tational time for the GRAPE system composed of multiple 
GRAPE 7, we need to use a limited rmmber of particles in 
each simulation. We run a model which requires ^ 100 CPU 
hours of the adopted GRAPE 7 system for at least ~ 6 Gyr 
evolution of a spiral in a group. If we use the same particle 
mass both for stellar particles in a spiral and for dark mat- 
ter particles of its host group, then we need at least ^ l(f 
particles for the dark matter. Numerically, it is very costly 
for the present work to run models with such a large particle 
number, given the limited CPU time. We thus use a large 
mass of ~ 10^ — lO^^M© for the dark matter particles of a 
group to dramatically reduce the total particle number of 
the group. 

Owing to this large mass, dynamical interaction be- 
tween group dark matter particles and a disk galaxy can 
cause unrealistic dynamical heating of the stellar disk (and 
possible transformation from spirals into SOs), which should 
be avoided in the present study. We therefore adopt a large 
gravitational softening length (e) for the dark matter par- 
ticles of a group to significantly suppress unrealistic tidal 
heating of a spiral by the dark matter particles: we use dif- 
ferent softening lengths for the different components in the 
simulation (e.g., dark matter halo and disk of a spiral and 
dark matter halo of its host group). We confirm that galaxies 
with large pericenter distances of their orbits with respect to 
their host group's center (thus with no/little tidal galaxy in- 
teraction) can not be transformed into SOs with thick disks: 
artificial tidal heating by group dark matter particles is well 
suppressed. 

The gravitational softening lengths for the dark mat- 
ter particles of a spiral, the stellar particles, and the dark 
matter particles of its host group are denoted as ea.d, ea.s, 



and Egr, respectively. We determine e for each component 
based on the mean particle-particle separation within the 
half-number radius of the spatial distribution of paxticle for 
each component. Furthermore, when two different compo- 
nents interact gravitationally, the mean softening length for 
the two components is applied for the gravitational calcula- 
tion. For example, e = (ea.s -|-egr)/2 is used for gravitational 
interaction between stellar particles in a disk and the dark 
matter particles of the host group. The gravitational soften- 
ing lengths for particles of the bulge, the gaseous disk, and 
new stars in a spiral are the same as ea.s, and e for group 
member galaxies (represented by point-mass particles) is as- 
sumed to be the same as egr. 

The total number of particles used in the fiducial model 
is 216,700 and 100,000 for a spiral and its host group, respec- 
tively, and ed,d, ed,s and £gr are set to be 1.79 kpc 0.25 kpc, 
and 10.78 kpc, respectively. As described above, these soft- 
ening lengths depend on i?d (i-c, smaller for less massive 
disks). The total number of dark matter particles in a group 
is linearly proportional to Mgr: it is 500,000 for the mod- 
els with Mgr = IO^^^Mq. We consider that as long as we 
investigate galaxies with Mgai ranging from ~ lO^^M© to 
~ IO^^Mq, dynamical friction of the galaxies against the 
dark matter halos of their host groups can be properly in- 
vestigated for such particle numbers used for groups. Also 
the total number of bulge particles is also linearly propor- 
tional to Mb (= /bMd) so that models with larger /b have 
larger bulge particle numbers. 

2. 5. 3 Main points of analysis 

We consider that if a disk galaxy has no spiral arms (by 
naked-eye inspection) at the final time step of each simula- 
tion, then the simulated galaxy can be classified as an SO 
galaxy. This is consistent with the canonical classification 
method that has been used in observational studies (e.g., 
Sandage 1961). We also use the smoothed two-dimensional 
(2D) stellar distributions of simulated galaxies to demon- 
strate more clearly that the simulated SOs have no spiral 
arms if they are projected onto the sky. The method to 
derived the smoothed 2D distributions is given in the Ap- 
pendix A. We investigate star formation histories, the mass 
fractions of stripped gas and stars, maximum rotational 
speeds of stars (Vm), velocity dispersions of gas and stars 
(cr), and evolution of the chemical abundances of the disks 
and nuclear regions. We also investigate the spatial distri- 
butions of the stripped gas and stars, which may well be 
identified as intragroup gas and stars. 

Our particular interest is in how spirals with different 
Hubble-types (i.e., different /b) evolve as a result of galaxy 
interactions in groups. The Hubble-types Sd, Sc, Sbc, Sb, 
and Sa (from late- to early-types in order) have /b = 0.02, 
0.10, 0.19, 0.32, and 0.70, respectively, in the present model 
(or B/T = 0.02, 0.09, 0.16, 0.2, and 0.41, respectively). 
For each model, we investigate the increase of the bulge 
mass (AMb) in a spiral due to galaxy interactions by esti- 
mating the total mass of new stars formed from gas within 
R ^ 0.2Jid. The growth rate of a spiral bulge is measured by 
AMb/Afb(0), where Mb(0) is the initial mass of the bulge. 
We do not adopt the same bulge-disk decomposition tech- 
nique as used for estimating B/T from the radial light pro- 
files of SOs in observations (e.g., Laurikainen et al. 2009). 




Figure 1. The time evolution of the stellar distribution of a spiral galaxy projected onto the x-y plane (left four panels) and onto the 
x-z plane (right four panels) at different times (T) for the fiducial model. The time T (in units of Gyr) is shown in the upper left corner of 
each panel. The solid blue line in the upper left panel on the left-hand side represents the orbit of the spiral for the last 5.6 Gyr. The red 
particles represent group member galaxies. The upper four panels on the left and right describe the time evolution of stellar distributions 
in smaller- and larger-scales, respectively. Stellar particles are not seen in the upper two of the left four panels, because tidal stripping 
can not happen before T = 2.8 Gyr. In order to clearly show the formation of intragroup stars, only stars that are stripped from the 
spiral (i.e., those located at /? > 2R^, where R is the distance from the center of the spiral) arc shown in the left four panels. 



This is mainly because we can not calculate the spectral 
energy distribution at each local region (for estimating the 
radial hght profiles) owing to the fact that the present nu- 
merical code has not been yet combined with the stellar 
population synthesis codes. In our future papers, we will 
combine our new chemodynamical code with a stellar pop- 
ulation synthesis one and thus be able to derive B /T in the 
exactly the same manner as observations: we will discuss 
this important issue extensively in a forthcoming paper. In 
the following, T in a simulation represents the time that has 
elapsed since the simulation started. 

3 RESULTS 

3.1 Fiducial model 

Figure 1 summarizes (i) the time evolution of the orbit and 
stellar distribution of a spiral in a group, and (ii) that of 
the distribution of group member galaxies, for the fiducial 
model. As the spiral approaches the pericenter of its orbit 
for the first time, the tidal field of the group starts to dy- 
namically heat up the stellar disk to some extent. The tidal 
field of the group, however, can only influence the outer part 
of the stellar disk of the spiral. In the early evolution of the 
spiral within the group (T < 2.8 Gyr), it only encounters 
low-mass dwarf galaxies so that galaxy interactions do not 
dramatically influence the evolution. However, as the spiral 
slowly sinks into the central region due to dynamical fric- 



tion, it more frequently interacts with group member galax- 
ies due to the higher number density of the galaxies their. 
This sinking of the spiral can not be seen in low-mass disk 
galaxies owing to much less effective dynamical friction for 
the spirals in the group. After T ~ 4 Gyr, the spiral experi- 
ences a number of tidal encounters with galaxies which have 
comparable masses (to the spiral) so that the stellar disk 
is strongly heated up to form a thick disk with a vertical 
velocity dispersion of ~ 60 km s^^ at i? ~ 10 kpc. 

Spiral arms within the disk can gradually disappear 
due to the strong tidal heating during the interactions with 
group member galaxies. About 24% (11%) of the old stars 
initially in the disk can be finally located at i? ^ 2i?d and 
\z\ ^ 0.1i?d [R < 2Ra and \z\ ^ 0.2iid). This stellar com- 
ponent can therefore be regarded as a stellar halo around 
the SO at T = 5.6 Gyr. The edge-on view of the disk galaxy 
becomes strikingly similar to those of SOs at T = 5.6 Gyr. 
About 11% of stars can be stripped to be located at i? > 2Rd 
during the morphological transformation of the spiral into 
an SO. These stripped stars may well be identified as intra- 
group stars drifting freely in the group. Thus the spiral is 
transformed into an SO with a thick disk and a stellar halo 
as a result of the multiple tidal interactions within the group 
within 5.6 Gyr. 

Figure 2 clearly shows that the simulated galaxy has 
a barred structure yet no clear spiral-arm structures of old 
stars in the simulated 2D image at T = 5.6 Gyr: this galaxy 
can be classified as an SO (or a barred SO). The new stars 
have a very compact distribution (see Figure Al in the Ap- 
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u (T=5.6 Gyr) 




Figure 2. The final smoothed 2D density distribution of old 
stars projected onto the x-y plane in the fiducial model. This 
logarithmic density map (/is) clearly demonstrates that the sim- 
ulated galaxy does not have spirals and thus can be identified as 
an SO. 




Figure 3. The time evolution of the star formation rate in the 
spiral for the fiducial model. 



pendix A), because star formation in the central region of 
the bulge is greatly enhanced in a repetitive manner dur- 
ing multiple tidal interactions with group member galaxies. 
About 60% of the initial gas can be rapidly consumed by star 
formation so that the final remnant can have a gas mass 
fraction of only ~ 0.04 (i.e., the formation of a gas-poor 
SO). The isolated model with no group tide and no galaxy 
interactions shows that only 20% of the initial gas can be 
converted into new stars after 5.6 Gyr of evolution. These 
results demonstrate that gas-poor SOs can be formed from 
multiple tidal interactions in groups. The final distribution 
of gas has a weak asymmetric feature (spiral-like morphol- 
ogy) , which is due to the last tidal interaction with a massive 
group member galaxy. 
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Figure 4. The time evolution of the bulge growth rate 
{AM\^/M^,{0)) for the fiducial model (dotted) and the isolated 
model (solid) in which neither group tidal field nor galaxy inter- 
actions are included. 
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Figure 5. Radial profiles of rotational velocities {V, solid) and 
velocity dispersions (cr, dotted) for old stars (red) and new stars 
(blue) in the simulated SO at T = 5.6 Gyr for the fiducial model. 
Here V and a are derived from line-of-sight velocities along the 
major axis of the simulated SO projected onto the x-z plane (i.e., 
edge-on). For this estimation, the jy-component of a velocity for 
each individual particle is used. 




Figure 6. Radial profiles of V (thick) and a (thin) at T = 
Gyr (blue solid), 2.8 Gyr (green dotted), 4.2 Gyr (magenta short- 
dashed), and 5.6 Gyr (red long-dashed) for the fiducial model. 
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Figure 7. The time evolution of tiie formation rate of new stars 
witliin tlie central 0.2R^ (i-e., bulge region) of the spiral (upper) 
and the mean metallicity for the new stars there (lower) for the 
fiducial model. 



Figure 3 shows that the spiral experiences a number 
of moderately strong "starbursts" (with the maximum star 
formation rate of ~ 6Mq yr~^) triggered mainly by tidal in- 
teractions with group member galaxies. The star formation 
rates in the intervals between the starbursts are also higher 
than the quiescent level before the spiral is influenced by 
the group tide (T < 1 Gyr), which means that interactions 
with the group member galaxies can also enhance the aver- 
age level of star formation. This enhancement is due largely 
to the significantly increased gas density in the inner region 
of the disk. We have run a model in which group member 
galax;ies are not included and other model parameter values 
are exactly the same as the fiducial model, and confirmed 
that such moderately strong repetitive starbursts as shown 
in the fiducial model do not happen in this model (also only 
35% of gas is consumed by T = 5.6 Gyr). This means that 
multiple tidal encounters are responsible for the significantly 
enhanced level of star formation in the spiral. 

Figure 4 shows that (i) the bulge of the spiral can in- 
crease more rapidly in its total mass compared to the iso- 
lated model, and (ii) the more rapid increase can be seen 
clearly at T > 2 Gyr when the spiral starts undergoing tidal 
interactions with group member galaxies. The final bulge 
mass is a factor of ~ 20% larger than its original value in 
the fiducial model. This rapid increase in the bulge mass 
results from the fact that a large amount of the ISM can be 
transferred to the bulge region in a repetitive mann(>r dur- 
ing multiple tidal interactions with group member galaxies. 
The bulge growth in the isolated model is due to the inner 
transfer of gas by dynamical action of a stellar bar in the 
disk (i.e., due to "secular evolution"). 

Figures 5 and 6 summarize the final stellar kinematics of 



the simulated SO and the time evolution of the kinematics 
during the morphological transformation. The kinematical 
properties are derived along the major axis of the stellar 
distribution of the SO projected onto the x-z plane and the 
y-component velocities of stars arc used for radial profiles of 
line-of-sight rotational velocities (V) and velocity dispersion 
(a). The rotational profile of the old stars (V) at T = 5.6 Gyr 
shows the maximum rotation (V^m) of 112km s~^, which is 
significantly smaller than the original Vm {~ 220kms~^). 
This is due largely to dynamical heating of the SO progenitor 
spiral by galgixy interactions with group member galax;ies. 
The velocity dispersion profile (cr) of the old stars is rather 
flat at |_R| > 4kpc and has a peak (cro) of 144kms~^ (thus 
Kn/cro = 0.78). 

Although the rotational profile of the new stars ap- 
pears to be more disturbed, there is no/little difference 
in the radial profile of V between old and new stars: Vm 
(= 128kms~^) of the new stars is slightly higher than that 
of the old stars. The new stars have a central dip in their 
a profile (i.e., ao = 92kms~^), which is seen in most mod- 
els in the present study. This smaller ao results from the 
fact that the new stars of the bulge are formed from gas 
after it loses a substantial amount of its kinetic energy due 
to gaseous dissipation (mostly by shock) during galaxy in- 
teraction. Given that the central new stars are young and 
metal- rich (as described later), these results mean that the 
younger, metal-rich populations in the SO bulge has differ- 
ent kinematics from the older, metal-poor populations in the 
bulge. 

The time evolution of the V and a profiles demonstrates 
that the spiral with a high rotation amplitude (Kn/co ~ 3) 
and a steep radial gradient of a can be transformed into 
an SO with a lower rotation amplitude (Kn/co ~ 1) and a 
fiat radial profile of a. It is confirmed that these kinemat- 
ical properties of the simulated SO can be seen in different 
projections (e.g., kinematics derived from Vz of stars in the 
x-y projection). Thus these results mean that if SOs arc the 
remnant of violent tidal interactions in groups, then kine- 
matical properties of SOs can be significantly different from 
spirals in the sense that SOs are dynamically "hotter" than 
spirals. 

Figure 7 shows that gas fueling to the central region 
{R ^ 0.2_Rd = 3.5 kpc) becomes more efficient after T ~ 
3.5 Gyr so that new stars in the central regions can be formed 
efficiently. The more efficient central star formation after 
T ^ 3.5 Gyr is due to more frequent tidal interactions with 
group member galaxies that are as massive as the spiral. Due 
to the infall of metal-poor gas from the outer part of the disk, 
the mean metallicity of the new stars can become slightly 
smaller at some time steps (e.g., T ~ 4.7 Gyr). However, 
such a slight decrease of the mean metallicity is a temporary 
phenomenon, and the metallicity can soon start to increase 
owing to chemical enrichment. Thus new stars in the bulge 
of the SO have an age-metallicity relation whereby younger 
stars are likely to have higher metallicities. 

3.2 Pcirameter dependences 

Although transformation processes from spirals into SOs 
are similar between the fiducial model and other models 
in which SOs are formed, the details of the transformation 
processes depend on model parameters. Also not all of the 
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Figure 8. The dependence of final morphological properties 
of the simulated SOs on the Hubble-types of the SO progeni- 
tor spirals. The results for the Sa, Sb, Sd, and LSB models are 
shown in upper left, upper right, lower left, and lower right, re- 
spectively. The models with Sa, Sb, and Sd Hubble-types have 
Md = 6 X IO^^Mq and = 17.5 kpc whereas the LSB model 
has Md = 6 X 1O^M0 and R^i = 10.0 kpc. 
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Figure 9. The dependence of (_B/T)j on (B/T)j for the simu- 
lated SOs in the representative 20 models with different /j, , /g and 
orbits in a group with Mgj. = 2 X 10^^ Mq. The dotted line rep- 
resents the border line where (B/T)^ = (B/T)-^ whereas the two 
horizontal dotted lines show the observed range of B /T that typ- 
ical SOs have. Here it should be stressed that the simulated B/T 
are not derived in the same way as done in recent observations 
(which uses bulge-disk decomposition methods). The increase in 
B/T in the present simulation is due to the mass increase of new 
stars in the simulated spiral galaxy. 



Figure 10. The distribution of stripped HI gas projected onto 
the x-y plane in a group with Mgr = 2 X 10^^ Mq for the low-mass 
LSB models with = 6 X IO^Mq, _Rd = 10.0 kpc, /t = 0.17, 
flini = 2rs, and /v = 0.25 (blue), with = 1.2 X IO^Mq, 
fld = 6.6 kpc, /b = 0.17, Rini = 4ra, and /v = 0.25 (red), with 
A/d = 1.2 X IO^Mq, Rd = 6.6 kpc, /b = 0.17, Rini = 6rs, and 
/v = 0.25 (green), and with Md = 1.2 X lO^M©, _Rd = 6.6 kpc, 
/b = 0.17, -Rini = 4rs, and /v = 0.5 (magenta). The gas particles 
located at _R > 2_Rd are regarded as being stripped from the 
spiral in each model in this figure. The big circles represent the 
final locations of the SO galaxies in the four models. 



present models with different model parameters show such 
morphological transformation into SOs. We illustrate here 
the derived dependences on the model parameters. 



3.2.1 Orbit 

Spirals in models with smaller /v (thus smaller pericenter 
distance rp and larger orbital eccentricity Bp) are more likely 
to be transformed into SOs for a given set of other parame- 
ters, because they can be more severely influenced by group 
tide and more frequently interact with group member galax- 
ies in the central regions of groups. Also spirals in models 
with smaller /v can experience stronger starbursts and thus 
grow their bulge components to a larger extent owing to a 
larger amount of gas transferred to the central regions and 
converted into new stars there. For example, the model in 
which /v = 0.75 and other model parameters are the same 
as those in the fiducial model, does not undergo the mor- 
phological transformation from a spiral into an SO seen in 
the fiducial model. Furthermore, models with smaller Rini 
are more likely to be transformed into SOs for a given set of 
other parameters. 
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3.2.2 Hubble-type 

In the models with A/gr = 2 x 10^^ and iiini = 4rs and 
/v = 0.25 and 0.5, spirals with different /b (thus different 
Hubble-types) can be transformed into SOs as a result of the 
disappearance of their initial spiral arms. However, bigger 
bulges in the models with larger /b can severely suppress 
the gas transfer to the central regions of disks so that sec- 
ondary starbursts are not so strong in the models. This is 
mainly because (i) stellar bars play a role in radial transfer 
of gas within disk, and (ii) the formation of strong bars dur- 
ing tidal interaction can be suppressed by the presence of big 
bulges. Figure 8 shows the edge-on views of SOs formed from 
spirals with different initial /b (i.e., Hubble-types). One of 
the common features of these simulated SOs are that they 
have thick disks and metal-rich stellar halos. 

3.2.3 Gas mass fraction 

Although the spiral to SO morphological transformation pro- 
cesses do not depend strongly on /g, star formation histo- 
ries are quite different between models with different /g for 
a given set of other model parameters. Due to higher star 
formation rates in the central regions of spirals, the bulges 
can grow more rapidly and to a larger extent in models 
with larger /g. For example, the model in which /b = 0.5 
and other model parameters are the same as those in the 
fiducial model shows {B/T)f = 0.22 (in comparison with 
{B/T){ = 0.18 for the fiducial model). These results imply 
that the bulge growth rates of spirals during morphological 
transformation into SOs in groups depend on the gas mass 
fraction of the spirals when they enter into group environ- 
ments. 

3.2.4 LSB or HSB 

In the present study, low-mass galaxies are assumed to have 
low mean surface stellar densities (LSBs) so that their evo- 
lution in groups can be different from the evolution of high- 
mass galaxies (thus HSB) described above. Irrespective of 
Mgr, spirals in models with smaller Mgai (~ 1.2 x lO-'^'^M©) 
are strongly influenced by both the group tide and by the 
galaxy interactions with group member galaxies. As a result 
of this, the outer parts of stellar disks in these less mas- 
sive spirals are more severely disrupted and a larger fraction 
(~ 30%) of the old stars arc tidally stripped to become in- 
tragroup stars. The remnants of this violent tidal interaction 
appear to be dwarf SOs (dSO). Figure 8 shows one example of 
this dSO formation for the model in which Ma = 1.2 x 10^ M© 
and other model parameters are the same as those in the 
fiducial model. Such dSOs can have compact stellar struc- 
tures in the central regions so that they can be classified as 
nucleated dSOs. Thus the present simulations show that less 
massive spirals are more likely to be transformed into SOs in 
groups. 

3.2.5 Group mass 

Spirals in more massive group models with Mgr = 5 x 
10^^ Mq are less strongly influenced by group tide and in- 
teracting galaxies in comparison with those in models with 
Mgr = 2 X 1O"M0 for the same Rini/rs (~ 4) and /v 



(0.25 — 1.0). As a result of this, spirals can continue their 
gradual star formation and keep their spiral-arm structures 
in the models: transformation from spirals into SOs due to 
tidal interactions is less likely in more massive groups for the 
above orbits. This is true for other massive group models 
with Mgr = IO^'^Mq. The main reason for this is that the 
relative velocities of two interacting galaxies are higher in 
more massive group models so that galaxy interaction can- 
not so strongly influence the evolution of spirals. It should be 
stressed, however, that if spirals are initially located in the 
central regions (< 2rs) in these more massive group models, 
then spirals can also be rapidly transformed into SOs due to 
the strong group tidal field and frequent galaxy interactions 
there. 

3.2.6 Inf all model 

In infall models, it depends on the impact parameter, by, 
and /v whether spirals are strongly influenced by group tide 
and galaxy interactions so that their evolution can be sig- 
nificantly changed in comparison to isolated evolution. For 
models with Mgr = 2 x 10^"^ M©, spirals cannot be influenced 
strongly by tidal effects in the group as long ashy > ra (even 
for /v = 0.2). Infalling spirals with by ~ 0.5rs and /v = 0.2 
can finally have very eccentric orbits within the host group 
so that they can be strongly influenced by the group tide 
and the group member gala^xies when they are in the cen- 
tral region of the group (i.e., at their pericenter passages). 
Due to the longer orbital periods of infalling spirals, the 
timescale for them to be transformed into SOs due to galaxy 
interactions is longer. 

3.2.7 Bulge growth rate 

Figure 9 shows the dependences of {B/T){ on {B/T)\ for 
the 20 simulated SOs in a group with Mgr = 2 x lO^M© 
for models with Md = 6 x lO^^M© yet different /b, /g, and 
orbits. Clearly the B /T can significantly increase in spirals 
with different Hubble-types (i.e., /b) during their morpho- 
logical transformation into SOs due to multiple galaxy in- 
teractions within the group. However, the final {B/T)f for 
spirals with (B/T)i < 0.2 cannot bo close to that typically 
observed (« 0.3 - 0.5; Christlein & Zabludoff 2004), which 
implies that typical SOs are not likely to be formed from 
later-type spirals (i.e., later than Sbc). Although early- type 
spirals can be transformed into SOs with bigger bulges, the 
original bulge sizes do not change during the transformation 
into SOs. 

3.2.8 Formation of intragroup gas 

Figure 10 summarizes the unique distributions of stripped 

gas from low-mass galaxies with Aid ranging from 1.2 x 
IO'^Mq to 6 x IO^Mq in a group with Mgr = 2 x IO^^Mq. 
These stripped gas distributions are well outside their ini- 
tial disks {R > 2J?d) and are seen to be drifting freely in 
the intragroup space so that they would be identified as 
intragroup HI gas. The structure and kinematics of the in- 
tragroup gas depend strongly on the orbits of their host 
spirals. Apparently isolated giant HI clumps with no/little 
optical counter parts can be seen in some models (e.g., in 



Formation of so galaxies 11 



the LSB model with Md = 6 x IO^Mq, Ca = 6.23, /b = 0.17, 
/g = 0.1, i?ini = 2ra, and /v = 0.25). The presence of these 
HI streams and isolated clouds in groups can be possible 
evidence for morphological transformation from spirals into 
SOs (as discussed later). 



4 DISCUSSION 

4.1 Formation of SOs with bulge growth 

One of the particularly intriguing results in recent observa- 
tions on SO evolution in groups and clusters of galaxies is 
that the number evolution of SOs is significantly more dra- 
matic in groups (or poor clusters) of galaxies with a < 750 
km s"^ at Q.l< z < 0.8 (Just et al. 2010). Other key obser- 
vational properties of SOs, which should be explained by any 
models of SO formation, is that the bulge fraction B /T (de- 
fined as the ratio of the bulge to total luminosity of a galaxy) 
is significantly larger in SOs than in spirals (e.g.. Dressier 
1980; Christlein & Zabludoff 2004). These two properties 
are not explained so simply by the gradual transformation 
scenarios of SO formation via ram pressure stripjiing (e.g., 
Abadi et al. 1999) and removal of galactic halo gas (e.g., Lar- 
son et al. 1980; Bekki et al. 2002; Bekki 2009), because such 
scenarios do not predict the significant increase of the mass 
concentrations of spiral galaxies. Previous numerical simu- 
lations of the interactions and merging between two galaxies 
showed that tidal interaction and merging can transfer gas 
from the outer parts of disks into the inner regions to grow 
the bulges (e.g., Noguchi 1988; Bekki 1998) and thus may 
well be the physical mechanism for SO formation. 

The present numerical simulations for evolution of spi- 
rals interacting with numerous group members have shown 
that multiple slow tidal encounters can be responsible for 
transformation from spirals into SOs with bulge growth. The 
growth rates of galactic bulges depend on the masses, initial 
B/T, and orbits of the original spirals for a given group, 
and thus not all of spirals can become SOs with big bulge 
{B/T > 0.5). Although spirals with Sd morphology-type 
can significantly grow their bulges owing to the formation 
of strong bars and subsequent gas infall onto the central re- 
gions (resulting from dynamical action of bar on ISM) dur- 
ing tidal interaction, the final B/T can be still smaller than 
the typical one for SOs. Thus SOs with large B/T (> 0.5) 
can hardly bo formed from later-type spirals with Sd, Sc. 
and Sbc Hubble-types via tidal interactions in groups. Other 
physical mechanisms that increase significantly bulge masses 
(e.g., minor and unequal- mass merging) are necessary for the 
formation of SOs with bigger bulges. 

The present study has shown that less massive spirals 
are more likely to be transformed into SOs, because they 
are more likely to interact with galaxies that are larger than 
themselves in groups and arc more susceptible to group tidal 
fields. This suggests that the observed rapid evolution of SOs 
in groups (e.g., Wilman et al. 2009; Just et al. 2010) is due to 
less luminous galaxies in groups. Also the present study has 
shown that spirals in the inner regions of groups can be more 
rapidly and dramatically transformed into SOs, which sug- 
gests that number evolution of SOs can be more clearly seen 
in the central regions of groups. Thus it would be impor- 
tant to investigate observationally how the number fraction 



of SOs depends on redshift in groups for a given mass-range 
and whether the dependences are different between the inner 
and outer regions of groups. 

4.2 Origin of the Tully-Fisher relation in SOs 

A number of authors have recently discussed the origin 

of SOs using the Tully-Fisher relation (TFR; Tully-Fisher 
1977) of SOs in different environments (e.g., Hinz et al. 2003; 
Bedregal et al. 2006; Williams et al. 2010). Hinz et al. (2003) 
found (i) a larger fraction of stellar mass in SOs than in late- 
type spirals and (ii) a larger scatter (a ^ 1 mag) of the 
TFR, in SOs than in spirals. They therefore suggested that 
SOs are likely to be formed from minor mergers rather than 
from simple truncation of star formation and subsequent 
disk-fading of late-type spirals. Bedregal et al. (2006) found 
that SOs lie systematically below the TFR for nearby spirals 
and have a larger scatter in the TFR than spirals and thus 
suggested that simple disk-fading alone cannot explain the 
observations so well. Williams et al. (2010) found that SOs 
are on average fainter in the Ks band at a given rotational 
velocity and suggested that the observed offset can not be 
simply explained by disk-fading of spiral galaxies. 

The present simulations have shown that if SOs are 
formed from spirals via tidal interaction in groups, then SOs 
can have smaller maximum rotational velocity (V^m) than 
their progenitor spirals owing to tidal heating of the disks. 
Although gas can be converted into new stars more effi- 
ciently during SO formation via tidal interactions, some frac- 
tion (10 — 30%) of old stars initially in the disks of spirals 
can be stripped to become intragroup stars. Therefore, it 
is unlikely that SOs have systematically larger disk masses 
in comparison with spirals. Our simulations show that the 
maximum circular velocities (K:,m) of SOs are very similar 
to or only slightly lower than those of their original spi- 
rals, though Vm of SOs are significantly lower than spirals: 
Vc,m = 245 km s^^ for the original spiral and 235 km s~^ 
for the final SO at T = 5.6 Gyr (due to mass loss) in the 
fiducial model. Our simulations therefore imply that SOs 
formed from tidal interactions can show Vm similar to or only 
slightly lower than spirals for a given luminosity a few Gyr 
after their formation (when the secondary starburst popula- 
tions formed during SO formation via tidal interaction have 
faded out). Observational results by Williams et al. (2010) 
do not show systematically lower K,m in SOs for a given lu- 
minosity, though their SOs are not within groups and their 
Vc,m values are derived from gas (not stars as done in the 
present simulations). 

As pointed out by many authors, there could be a num- 
ber of viable physical mechanisms for SO formation (see §1 
in this paper). If SOs are formed from disk- fading after trun- 
cation of star formation due to removal of galactic halo gas 
(e.g., Larson et al. 1980; Bekki et al. 2002), then SOs are 
likely to lie systematically below the TFR of spirals (e.g., 
Bedregal et al. 2006) . If SOs arc formed from spirals by tidal 
interaction with some fractions of original dark matter and 
disk stars being stripped, then SOs can show slightly lower 
Vc,m for a given luminosity or mass owing to the mass loss 
during tidal interaction. Therefore, the observed large scat- 
ter in the TFR of SOs (e.g., Hinz et al. 2003; Williams et 
al. 2010) appears to be consistent with SO formation via 
different physical mechanisms. It is currently impossible to 
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discuss whether the TFR of SOs in groups, where galaxy 
interactions can be one of the major mechanisms for SO for- 
mation, is different to that in the field and dusters due to the 
lack of observational data sets. Thus future observations on 
the TFR in different environments will enable us to discuss 
what mechanism(s) dominates SO formation in the field, and 
in groups and clusters, based on the locations of SOs in the 
TFR. 

4.3 Formation of intragroup cold gas 

Recent observations have confirmed that nearby clusters 
have intraclustcr planetary nebulae (e.g., Arnaboldi ct al. 
2003; Arnaboldi 2010 for a recent review) and globular clus- 
ters (e.g. Bassino et al. 2003; Lee et al. 2010), which implies 
that these intraclustcr stellar objects originate from tidal 
stripping of stars and globular clusters from cluster mem- 
ber galaxies. Numerical simulations based on ACDM models 
showed that such intraclustcr staxs can be formed during the 
hierarchical assembly of clusters of galaxies and the simu- 
lated radial distributions of the stars are consistent with 
the observed ones (e.g., Murantc ct al. 2007). Furthermore, 
cosmological N-body simulations with a formation model 
of globular clusters demonstrated that intraclustcr globular 
clusters can be formed from tidal stripping of globular clus- 
ters initially in low-mass galaxy-scale halos (Yaliagi & Bekki 
2005; Bekki & Yahagi 2006). Although these observational 
and theoretical studies have discussed extensively the origin 
of intraclustcr stellar objects in clusters, they did not discuss 
the possible existence of intragroup gas. 

The present study has demonstrated that if SOs arc 
formed from gas-rich spirals via tidal interactions in groups, 
then there should be intragroup HI gas that can have unique 
distributions like very long tails, broken-rings, and appar- 
ently isolated massive clouds. Therefore it would be reason- 
able to claim that the presence of such intragroup HI gas 
is possible evidence for morphological transformation from 
spirals into SOs via tidal interactions. Also if SOs are forming 
in groups at < 2: < 0.8, as suggested by recent observations 
(e.g., Just ct al. 2010), then a larger amount of intragroup 
HI gas should be observed in groups at such redshifts. The 
present simulations also predict that tidal stripping of gas 
is more efficient in the outer parts of disks (where gaseous 
metallicities arc lower owing to the negative metallicity gra- 
dients) for less massive galaxies so that the intragroup gas 
can be composed mostly of metal-poor gas. This mctal-poor 
gas would finally mix with the hot gaseous halos of groups 
and consequently change the metallicity of the halo gas. 



5 CONCLUSIONS 

We have numerically investigated in detail how spirals dy- 
namically interact with group member galaxies and group 
tidal fields and consequently evolve into SOs using chemo- 
dynamical simulations of galaxies. We have analyzed struc- 
tural, kinematical, and chemical properties of the simulated 
SOs in groups for a variety of model parameters (e.g., Md, fh, 
and Afgr) that can the control physical properties of galaxies 
(e.g., Hubble- types) and groups. We summarize our princi- 
ple results as follows. 

(1) Spirals can be transformed into gas-poor SOs with 



no remarkable spiral arms within ~ 6 Gyr evolution of the 
spirals within groups with Mgr = 2 x lO^'^Mo. Multiple slow 
tidal encounters with group member galaxies can dynami- 
cally heat up the disks of spirals so that their initially thin 
disks can be transformed into thick ones during morpholog- 
ical transformation from spirals into SOs. Such tidal interac- 
tion can also trigger moderately strong (up to ~ IOMq yr~^ 
for /g = 0.1) starbursts in the inner regions of their bulges so 
that bulges can grow significantly during SOs formation. The 
strong tidal fields of groups can also play a role in morpho- 
logical transformation from spirals into SOs for less massive 
spirals. 

(2) The details of the processes responsible for the mor- 
phological transformation of spirals into SOs depend strongly 
on galaxy masses, the orbits of spirals within groups, and 
the total masses of groups (Mgr). Spirals initially located in 
the inner regions of groups can be more rapidly converted 
into SOs due to the larger galaxy densities and stronger tidal 
fields of groups there. Less massive spirals are more likely 
to be transformed into SOs, because they can interact with 
galaxies that are significantly more massive than themselves 
and are susceptible to tidal fields of groups. 

(3) Due to significantly enhanced star formation and ef- 
ficient chemical enrichment in the inner regions of spirals 
during morphological transformation into SOs, the bulges of 
the final SOs have young and metal-rich stellar populations. 
The young, metal-rich stellar populations are located in the 
inner regions of SO bulges and have dynamically colder (i.e., 
smaller a) than old stellar populations initially in the bulges. 
The mass fractions of younger and metal-rich populations in 
SO bulges depends strongly on the initial gas mass fractions 
(/g). Less massive spirals (Aid ~ 1O""'M0) with small bulges 
can be transformed into nucleated SOs due to strong nuclear 
starbursts occurring during the morphological transforma- 
tion process. 

(4) Stellar disks of SOs formed from tidal interactions 
in groups have significantly fiatter radial profiles of line-of- 
sight velocity dispersions (cr) in comparison with original 
spirals, because the outer parts of the original disks (where 
a are initially small) can be dynamically heated up more 
strongly. The maximum rotational velocities (Kn, which is 
significantly lower than the maximum circular velocity, Vc,m) 
of disks can become significantly smaller during morphologi- 
cal transformation owing to dynamical heating by slow tidal 
encounters and group tidal fields. Thus, if SOs are formed 
from tidal interaction in groups, then they can have smaller 
Vm for a given luminosity. 

(5) Intragroup cold HI gas can be formed from gas 
stripped from spirals being transformed into SOs due to tidal 
interactions in groups. Such intragroup gas can have unique 
spatial distributions (e.g., rings and very long tails) that 
reflect the orbits of gas-rich spirals from which the gas orig- 
inates from. Some intragroup gas clouds are far away from 
the galaxies they originated from so that they can be identi- 
fied as apparently isolated massive HI clouds with no optical 
counterparts in groups. Since most intragroup gas is from 
the outer parts of disks in less massive spirals, the metallic- 
ity of the gas should be significantly lower. The presence of 
intragroup HI gas with unique spatial distributions is sug- 
gested to be possible evidence for morphological transforma- 
tion from spirals into SOs due to tidal interaction in groups. 

Finally, tidal interactions with group member galaxies 
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can be just one of a number of possible physical mecha- 
nisms for the transformation of spirals into SOs. SOs with 
thin disks can be formed from the gradual truncation of star 
formation caused by halo gas stripping whereas SOs with big 
bulges {B/T > 0.5) and with thick disks can be formed from 
minor and unequal-mass galaxy merging. In our forthcom- 
ing papers, we discuss the physical properties of SOs formed 
by these physical processes that were not explored in the 
present paper. 
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APPENDIX A: DERIVATION OF 2D DENSITY 
FIELDS 

In order to compare the smoothed 2D density fields of simu- 
lated SOs with the observed images of galaxies (on the sky) 
in a more self-consistent manner, we adopt the smoothing 
method that is identical to that used for analysis of the simu- 
lated stellar distributions in the remnants of galaxy mergers 
(Bekki & Peng 2006). We therefore briefiy summarized the 
adopted method in the present paper. At the position of 
each stellar particle, we apply a local linear smoother us- 
ing a Gaussian kernel function with the smoothing length 
of 0.05 in our units (corresponding to 0.875 kpc). We divide 
the relevant region with the size of ~ 30 kpc (as shown in 
Fig. 2) into 50 x 50 cells for a model in each projection 
and estimate the local smoothed stellar density (Es). We es- 
timate the 2D logarithmic density map {fis ~ log^QEs) for 
galaxy in a model in order to mimic the observed image of 
the galaxy on the sky. Total number of cells in each pro- 
jection is fixed at 2500 for all models in the present study, 
because this number is enough to show clearly the presence 
or absence of the spiral arms in the simulated disk galaxies. 

Figure Al shows the spatial distributions of old stars, 
gas, and new stars at T = 5.6 Gyr in the fiducial model, 
which are used for creating smoothed 2D stellar distribu- 
tions shown in Figure 2. Although an arm-like gaseous struc- 
ture can be discernibly seen in Figure Al, stellar distribu- 
tions are overall featureless so that the simulated galaxy 



with no spiral arms can be classified as an SO. Figure A2, 
which shows the smoothed 2D density profile for new stars, 
demonstrates that (i) the young component (formed from 
secondary starbursts) of the SO has a different spatial dis- 
tribution in comparison with the old one and (ii) it also 
has no clear spiral arm-like structures in the x-y projection. 
The central young component might well be observed as 
a "compact blue bulge" in observations. Figure A3 shows 
the smoothed stellar distribution of old stars viewed from 
edge-on (i.e., x-z) and confirms that the simulated galaxy 
appears to be "lenticular" like observed SOs. It is clear from 
this figure that the SO has an extended stellar halo that is 
formed from tidal stripping of stars during tidal interaction 
in the group. In order to clearly show how spirals are trans- 
formed into SOs due to tidal interaction, we have produced 
2D images on the time evolution of smoothed 2D stellar dis- 
tributions of galaxies: Figure A4 shows an example for this 
in the fiducial model. 



APPENDIX B: BULGE-DISK 
DECOMPOSITION 

Figure Bl shows the radial mass-density profiles separately 
for bulge and old stellar disk at T = in the original spiral 
galaxy for the fiducial model. This is simply a radial den- 
sity profile and not a luminosity-density profile: the radial 
luminosity-density profile depends strongly on the age and 
metallicity distributions of stars in the bulge and disk com- 
ponents. Figure B2 shows the final radial mass-density pro- 
files separately for the bulge, old stars, and new ones of the 
simulated SO at T = 5.6 Gyr in the fiducial model. The over- 
all profile for the total mass-density is significantly changed 
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Figure A3. The same as Figure 2 but for the x-z projection. 
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Figure Bl. Tlie final radial mass-density distributions (Sa(-R)) 
of bulge stars (solid blue), old disk ones (dotted red), and total 
(short-dashed green) projected onto the x-y plane at T = Gyr in 
the fiducial model. Here describes the initial disk size (=17.5 
kpc) in the fiducial model. 
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Figure B2. The final radial mass-density distributions (Ss(-R)) 
of bulge stars (solid blue), old disk ones (dotted red), new stars 
(short-dashed green), and total (long-dashed magenta) projected 
onto the x-y plane at T = 5.6 Gyr in the fiducial model. 



the origin of SOs in the context of the observed structural 
differences of stars between spirals and SOs in a fully self- 
consistent manner 



Figure A4. The same as Figure 2 but for different four time 
steps, T = 1.4 Gyr, 2.8 Gyr, 4.2 Gyr, and 5.6 Gyr. For con- 
venience, the color bar shown in Figure 2 to indicate physical 
mass-densities is not shown in each frame. 



between T = Gyr and T — 5.6 Gyr owing to tidal strip- 
ping of stars and secondary central starbursts during mor- 
phological transformation, though the central density does 
not change so much. The stellar bar formed in the stellar 
disk during tidal interaction can interact dynamically with 
the original bulge so that the radial density distributions of 
the bulge and disk components can change during the 5.6 
Gyr dynamical interaction in the group. In our future pa- 
pers, we try to produce radial hght profiles of the simulated 
SOs by combining the present chemodynamical code with 
the stellar population synthesis one so that we can discuss 



